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1. Introduction 

Unlike conventional, low-temperature superconductors high-Tc superconductors ex- 
hibit strong fluctuations due to the combined effect of the small coherence length ^, 
anisotropic layered structure and high temperatures. It has been argued that the 
Abrikosov vortex lattice melts as a consequence of these enhanced thermal fluctuations 
Ij^-]^. Quite recently there has been experimental evidence [p|-p!0| that in clean crystal 
samples (in the absence of twin boundary pinning) the Abrikosov flux lattice melts via a 
first-order phase transition In the flixx liquid state the vortex lines are free to move 
through the sample (except for their mutual repulsion) and will collectively drift in the 
presence of an external transverse current. This flux line motion will then in turn generate 
a finite voltage and lead to a nonzero linear resistivity. 

The phase diagram is changed if one includes the effects of disorder. Depending 
on the type and strength of the disorder, the vortex liquid state will persist down to an 
irreversibility line associated with a possible second-order phase transition to a vortex glass 
|T^ , or Bose glass state The melting of the flux lattice has clearly very important 



consequences as most vividly illustrated in fig. |T] where the mean-field phase diagram 
is contrasted with the phase diagram which includes effects of thermal fiuctuations and 
disorder. 

Because the flux- line liquid phase occupies a large portion of the H — T phase diagram, 
much of the efforts have been directed toward a better understanding of the properties of 
this phase in the presence of disorder. There has been considerable progress in under- 
standing the static properties of flux lines in the liquid phase The resulting phase is 
well described by a collection of directed flexible lines with a line tension related to Hd- 
The lines traverse the sample in the direction of the applied magnetic fleld and at flnite 
temperatures wander throughout the sample analogously to the Brownian motion executed 
by atoms or small molecules in conventional isotropic liquids. This flux-line liquid in the 
presence of point disorder is depicted in flg. |2| The flux-line states depicted in flg. have 
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been extensively studied in many experiments. Early experiments used the Bitter tech- 
nique in which the location of the flux line ends emerging from the sample is resolved by 
sprinkling magnetic powder on the sample surface. P,p^ The disappearance of the reg- 



ular, hexagonal pattern as the fleld or temperature are increased is suggestive of melting 
of the Abrikosov flux lattice. These experiments only directly probe the surface conflgu- 
ration of flux lines and therefore do not exclude a possibility of melting conflned to the 
surface. Indeed theoretical analysis shows that the surface interaction always dominates 
in determining the decay of translational correlations in the asymptotic long-wavelength 



limit. However, such large length scales have not been probed by the decoration 

experiments. Later vibrating reed experiments by Gammel et al.[^ have found a signal 
suggestive of the melting transition. Although very difficult, more direct measurements of 
the bulk properties of the flux-line liquid are possible. In principal the structure function of 
the interacting line liquid can be measured using neutron scattering techniques. These ex- 
periments can directly probe the density correlations in the line liquid, which are governed 
by very different physics than ordinary liquids of point particles with an isotropic structure 



function. fTq] Also, recently, new revolutionary electron holographic techniques have been 
used to image the motion of flux lines in thin Pb fllms. This probe, which can image flux 
lines in real time, can be used to study mixed states also in high-Tc superconductors, and 
can provide envaluable information about the dynamics of flux lines, fl^ 

The interactions between lines, thermal fluctuations and the effects of disorder are 
clearly very important and are the main physical features that must be taken into account 
by the theory. Signiflcant progress has been made by mapping the problem of flux lines 
onto the quantum statistical mechanics of interacting 2iI>-bosons, where the roles of h, 
temperature and mass of the bosons are played by the temperature, the inverse sample 
thickness and line tension, respectively. |^ With this mapping much of the insight gained 
from the study of systems like helium was taken over to the problem of line liquids. When 
the lines are long, (i.e. the sample along the applied fleld direction is thick), the temper- 
ature is high, and the 2L>-density is large, such that a typical line wandering distance is 



larger than the average interhne spacing, a highly entangled line liquid results. However, 
the entanglement should persist only if the line crossing barriers are significantly larger 
than the thermal energy. Although to date no detailed analysis of line crossing barriers 
exists, simple estimates give Ux/kBTmeit ~ 10 — 30, which translates into very slow relax- 
ational rates. The entangled phase corresponds to the Bose-condensed phase in the boson 
picture. 

The static structure function for the interacting fiux lines has been previously com- 
puted within the Bogoliubov approximation taking advantage of the boson mapping 
The contours of constant scattering intensity form a butterfiy pattern with two peaks, and 
is quite different from the structure function of isotropic fluids of particles. At long wave- 
lengths, away from any critical transitions, the theory of the vortex liquid can very well be 
described by a "hydro dynamic" model in terms of density fields, with phenomenological 
nonlocal coefficients. The advantage of starting with a microscopic description, however, 
is that the long wavelength description of the liquid can be understood in terms of a more 
basic, microscopic model, thereby providing a more detailed understanding. 

Although the boson mapping has been instrumental for understanding the static long 
wavelength behavior of line-liquids, unfortunately, there does not appear to be an obvious 
extension of this mapping to study the real hydrodynamics. Some progress has been 
made through phenomenological approaches in which the dynamic equations of motion 
are written directly for the coarse-grained density fields, using the static free energy, with 
phenomenological nonlocal coefficients determined by the static structure function. 

As for the statics it is useful to obtain the description of the dynamics, starting with 
equations of motion for individual interacting fiux lines, and to derive the dynamics for 
the observable hydrodynamic quantities like density. The goal is to construct a kinetic 
theory of flux-line hydrodynamics analogous to that of point liquids, which was useful in 
understanding hydrodynamics of simple liquids many years ago. With this approach it 
should be possible to calculate the hydrodynamic parameters (like line liquid viscosity, if 
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it exists) which arise completely from the fliix line interaction and entanglement effects 
and from the single line microscopic friction. The bare diffusion parameter will be an input 
to the theory, and is related to the real microscopic coupling of flux lines to the underlying 
crystal lattice, the Bardeen-Stephen friction coefficient. 



A question of renormalization of the bare diffusion coefficient of a tagged line by 
the presence of the flux-line liquid and the interaction with this liquid through excluded 
volume interaction is also of interest and is related to the line fluid viscosity. It is expected 
that if disorder is strong enough and the system is in an entangled regime a localization 
phenomena will take place, driving the renormalized diffusion constant to zero. 

In this paper we take the first step toward a description of the flux-line liquid in terms 
of a kinetic theory of line liquids. We introduce a formalism that is useful for microscopic 
calculations of the dynamics in the flux-line liquid phase. The hydrodynamics of the flux- 
line liquid is studied by starting with the microscopic description of the interacting flixx 
lines in terms of the Hamiltonian that includes the repulsive interline interactions and in 
the presence of quenched pinning disorder that couples to the density of lines. We expect 
that entanglement effects are in principle automatically incorporated in the full theory 
derived with this kinetic approach. The repulsive interaction will inhibit the lines from 
passing through each other and for high line densities will result in slowing down of their 
dynamics due to these constraints. It is not clear, however, what simplest approximation 
to the resulting interacting field theory will retain these effects. 

This paper is organized as follows. In Sec. 2 we introduce our microscopic model for 
the dynamics of interacting lines in the presence of quenched disorder, and in Sec. 3 for- 
mulate the dynamics in terms of a more convenient Martin-Siggia-Rose (MSR) description. 
In Sec. 4 the hydrodynamic description is discussed and the microscopic model is partially 
recast in terms of the density fields. In Sec. 5, by using the method of auxiliary fields, we 
integrate out the microscopic degrees of freedom and derive the effective MSR dynamic 
functional thereby obtaining a hydrodynamic description of the interacting flux-line liq- 
uid. In Sec. 6 we approximate this theory by truncating the expansion of the hydrodynamic 



functional at quadratic order, and within this approximation calculate the interacting dy- 
namic structure function in the presence of disorder. We analyze this dynamic structure 
function in various regimes in Sec. 7, and derive the corresponding static structure function 
demonstrating that it agrees with the result obtained via the boson mapping method. We 
find that on time scales longer than Rouse time (time required for the single line excitation 
of size L to relax elastically) or equivalently on length scales larger than the transverse 
line wandering length, the noninteracting (single fiux line) dynamics is dominated by the 
center of mass mode with a /c^ relaxational rate. In this regime we find that the interacting 
structure function reduces to that of hydrodynamic frictional diffusion, consistent with the 
phenomenological model of Marchetti and Nelson. |T^ The interactions between the lines 
additively renormalize the relaxational rate generating a crossover between the noninter- 
acting and interacting dynamics. This crossover occurs at fiux line length L = L/, and is 
physically related to the entanglement length defined in Ref. |l] . On wavelengths smaller 
than the transverse wandering length and for times shorter than the Rouse time we find 
that the noninteracting dynamics is controlled by the internal modes. In this regime we 
obtain a complicated wavevector-dependent renormalization of the dynamics summarized 
by the interacting structure function. In Sec. 8 we take the phenomenological approach to 
the hydrodynamic description of the fiux-line liquids and compare with the results of the 
kinetic approach derived in Sec. 7. Appendix A describes an independent derivation of the 
static structure function using the methods of auxiliary random fields. In Appendices B 
and C we derive the nonlinear terms in the expansion of the hydrodynamic functional and 
"hydrodynamic" Hamiltonian, and analyze the single line dynamics in various regimes. 



2. Dynamical Model 

The statistical mechanics of fiux-line liquids is very different from that of a liquid of 
point vortices because the lines are long and connected. Compared to these important 
topological properties the detailed internal structure of an individual fiux line is relatively 
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unimportant. The essential physics of the flux-hne hquid can therefore be described by the 
conformation and position of each hne. These configurations of vortex hues are charac- 
terized by a set of N functions Ili{z) = {ri{z), z), where ri{z) specifies the position of the 
ith hne in the (x, ?/)-plane as it wanders afong the direction of the apphed magnetic field 
H\\z through the sample of thickness L (see Fig. 2.1). The probability of an equilibrium 
configuration of interacting lines in the presence of disorder is given by the Boltzmann 
weight exp(— 7Y//cbT) with 

N L N L N L 

^=9^/ dz{d,nf + -J2 / dzV{n{z)-r-{z))+J2 dzU{n{z),z). (2.1) 

The first term describes the elastic energy of non-interacting lines with the line tension 
e, which for isotropic superconductors is given in terms of the London penetration length A 
and the ratio k = X/^ by e = (^o/47rA)^ In with the fiux quantum (pQ = hc/2e = 2 x 10~^ 
G cm^ and ^ the superconducting coherence length and the vortex-line core thickness. Here 
we are working in the regime for which e can be approximated by a constant, although 
with our formalism we can easily treat the case of a nonlocal elastic energy. 

The anisotropic superconductors can be well described by an effective mass tensor 
diagonal in the coordinate system with the z-axis aligned with the c-axis of the crystal, 

/Ml \ 

Mr^m = Ml . (2.2) 

V M3/ 

For these anisotropic layered compounds the line tension is considerably smaller e = 

where Mi is the in-plane anisotropic mass, and M3 ^ lO^Mi is the much larger 
effective mass describing the weak Josephson coupling between the planes. The above 
formula for e applies when the fiixx lines are dense (noA^ ^ 1, where A^ is the in-plane 
London penetration depth, and uq the average density of vortex lines per unit area). In 
the opposite limit noA^ ^ 1 the electromagnetic coupling between the planes is important 



and one gets e = e/ In k [p2| , p!2 



The second term in Eq.(2J.) incorporates the flux hne interactions. We treat the 
regime in which the hne coordinates vary slowly with although with our formalism we 
can easily extend our treatment beyond this regime. This leads to the interaction energy 
which can be expressed in terms of a pair potential which is local in z, and in the London 
limit is given by 

= ^ [^o(|r1/A) - Ko{\f\/0] , (2.3) 
where Kq{x) is the modified Bessel function with the asymptotics, 

In Eq.( p.3|) we have introduced a short distance cutoff, the superconducting coherence 
length ^. 

The final term in Eq. ( |2.1| ) is the contribution of the pinning impurities to the free en- 
ergy of the vortex line liquid. For simplicity we will take the quenched disorder interaction 
U{r) to be Gaussian with zero mean, 



t/(f, z) = , (2.5a) 



U{fi,Zi)U{f2, Z2) = F{fi - f2, Zi - Z2) , (2.56) 

where F{fx — r2^z\ — zq) encodes the strength and range of disorder correlations. We 
will later specialize to point, line and plane disorder. These three cases appear to be the 
most relevant experimentally as we discuss below, with oxygen vacancies and interstitials 
playing the role of point disorder, columnar defects and grain/twin boundaries as the line 
and plane disorders. 

The model introduced above has been used, in somewhat more specialized form to 
describe the static features of the vortex liquid. Nelson et al. used the boson mapping 
to compute the interacting static structure function in the dense phase where the density 
fiuctuations are small and mean field theory is a good description. Using renormalization 
group methods, they were also able to treat the dilute line liquid near iJci, where the 



fluctuations are strong. By matching to the dense phase, where the mean-field theory 
is accurate they computed the fiuct nation- corrected constitutive relation, B{H), and the 
static structure function near Hd- This work was further extended to include the effects 
of point pinning disorder on the static structure function in this phase. ||2^ 



Since we are seeking a hydrodynamic description of the fiux line dynamics we rewrite 
the interaction and disorder terms in a convenient form, 

Jo Jr„r„z 

V{fi - f2)n(fi, z, t)n(f2, z, t) - NLV{0) , (2.6a) 



' ri,r2,z 

^ dz U{n{z,t)) = / d^rdz U{f)J2^^^Hr-riiz,t)) , 
i=i J^ J i=i 

= J d'^rdzU{r}n{f,z,t) . (2.66) 

In above we defined /^^ = / d'^rdz. The self-energy term NLV{0) in Eq.( |2.6a| ), appropri- 
ately cutoff by ^, can be absorbed into the line tension energy, and we will therefore ignore 
it in the following analysis. The form of Eqs.(2.6 ) suggests that interactions in the liquid 
phase are naturally described in terms of the line density 7i(r, 2, t) = Xli^i ^^'^^ ^j(^5 ^))- 
Here we use the model defined by Eqs.( p.l|) , (2.5 ) to study the dynamics of the fiux- 
line liquid. We take the kinetic theory approach and write down the microscopic dynamic 
equations for each of the interacting lines. We assume that at this basic level the line 
fiuctuations are over damped and are therefore governed simply by Model A type Langevin 



dynamics ||2J] 



dtf,{z,t) = -Dj^^ + Uz.t) . (2.7) 
For simplicity we assume that the noise is Gaussian with zero average, and covariance 

{Q{z,t)C''{z\t')) = 2DkBT5,,5ab5{t-t')5{z-z') , a, 6 =1,2. (2.8) 
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The parameter D is the microscopic kinetic coefficient proportional to the inverse of 
the Bardeen-Stephen friction coefficient, p!9f 



70 = ^^'^^ ' (2.9) 

with as the superconducting coherence length in the copper-oxide planes and an is 
the normal-state conductivity. D describes the effective drag on a flux line due to the 
interactions of the normal core electrons with the underlying solid. In the simplest case 
of center-of-mass-dominated diffusion of noninteracting flux lines D is proportional to 
the macroscopic diffusion coefficient (see Sec. 7). The above dynamics of flux lines is quite 
similar to the dynamics of polymer melts, with the important difference that for flux lines 
there is no solvent. Since the vortex lines are not moving in any fluid solution. Rouse rather 
than Zimm dynamics applies, and D is a constant. However, we expect significant 
changes in the hydrodynamic parameters of the line liquid coming from the pinning disorder 
and flux-line entanglement, when there are significant barriers to line crossing. 

In Eqn. ( P77| ) we have for simplicity ignored the component of the line drag normal to 
the fiux velocity v, which can be accounted for by adding Az x dtri{z, t) to the left hand 
side of Eqn.( p77D . This term determines the Hall angle Oh, according to tan6'H = A.[p!9| 
and is generally quite small. 



3. Martin-Siggia-Rose Formulation of the Model 

We would like to derive a hydrodynamic description of the line liquid in terms of 
coarse-grained variables like the density of the flux lines, starting with the microscopic 
model presented in Sec. 2. This procedure is most conveniently implemented using Martin- 
Siggia-Rose (MSR) formalism [ p6| , p7| which allows the solution to the Langevin equation 



to be formulated in terms of a constrained path integral. 
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The idea is that instead of solving the Langevin equation for the conformation variables 
ri{z, t) in terms of the random forces C,i{z, t) and then computing the correlations functions 
by averaging over the noise with the Gaussian weight 

1 



^[Ci{z,t)] ~ exp 



dt / dz\Uz,t)Y 



(3.1) 



ADUbT 

one can consider fi{z,t) as the basic stochastic field with a path probability density 
W[ri{z,t)] defined by 



W[n{z, t)]Vn{z, t) = u[Uz, t)]VQ{z, t) 



(3.2) 



and eliminate the random forces in favor of the conformation variables. This is accom- 
plished via a constrained path integral over the noise with the Langevin equation as the 



constraint. [P6|,|28|,P9| 



To implement the procedure of MSR we note that the noise average of any observable 
C[ri(^,t)], with flux line conformational variables ri(z,t), as the solution of the Eq.( |2.7|) , 
can be expressed in terms of a constrained path integral, 



mm 



VUz,t)Vf,{z,t)J[fi\ Yl 5 

a,i,z,t 

1 



X 0{fi) exp 



dzdt\Qiz,t)\' 



(3.3) 



ADkeT 

Here / T>Q{z,t)Vfi{z,t) denotes a path integral over the noise and the conformation of 
the flux lines with the implied discretization of z and t to define the path integral. The 
quantity Jfr^] is the Jacobian of the transformation from Ci(-2, t) to Viiz, t) imposed by the 
functional 5-function. It ensures that the path probability density, W[fi{z, t)], with which 
the averages are computed is still normalized to 1, i.e. (1) = 1. 

We eliminate the functional 5-function by performing the integral over Q{z,t) and 
obtain 



{Oin)) = / Vn{z,t)J[n] 



X 0{fi) exp 



(3.4) 



11 



Further, it is convenient to perform a Gaussian transformation in order to "linearize" the 
argument in the exponential, the dynamic functional. This is accomplished by introducing 
auxiliary fields ri{z,t), usually called the response fields due their utility in computation 
of response functions as will become clear below, (see Appendix C) 



(3.5) 



where J^id[ri, ri] is the dynamic functional for a particular realization of disorder. 



dt I dz 





ftiz,t)DkBTrtiz,t) 



(3.6) 



It is convenient to work in Fourier representation. We reserve q for the wavevector in the 
2;-direction, k for the transverse wavevector and w for the frequency variable (see below). 
Since the fiux line ends are freely fiuctuating, these boundary conditions are naturally 
satisfied by the discrete cosine- Fourier transform with fi[z,t) and f'i{z,t) given by. 



"o{uj) + 2 ^ cos{qnz) r{qr. 

9n>0 



e^'"* cos 



((/„^) 2(^-0-1) r(Q„, a;) . 



(3.7) 



where q^ = rni/L and we have been careful to separately treat the center of mass, q^ ~ 
0, mode. This seperation will turn out to be essential in order to recover the correct 
hydrodynamic result at long wavelengths (see Appendix C). 

To examine the explicit form of J7id[ri,ri], we split it up into three contributions 
corresponding to the three terms of the Hamiltonian in Eq. (|2.1|) , 

Jid[ri-: Ti] = Jo[fi, Ti] + Ji[fi, fi] + Jd[fi, Ti] . (3.8) 

(i) iTof'^ij^i] is the dynamic functional for the noninteracting flux lines, 

Jo[n. r,] = L [ Yl [^^'(^-' ^)DkBTfnQn. -^) + r^Qn, c^) + iDeql) r^qn. -co)] 2'"'"-' 



'9"''^ a,i 



Ri{qn,Uj) ■ G{qn:Uj) ■ Rfiqn: "W) 



(3.9) 
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where k oj ^ Xln / '^^^/(^^r)^/ duj/{2Tr), and we defined, 



Ri{qn,uj) = (ri{qr,,uj),r'i{qn,uj)^ , (3.10a) 

(ii) J'i[ri, Ti] is the contribution to the total dynamic functional due to the interaction 
between fliix lines, 



N r r 

jS.r^] = -DTT dt dz {k-f,)V{k)e'''<^^-'^^^ , (3.11) 



where V{k) is the Fourier transformed flux-line interaction, Eg. (2.3) 



V{k) = / d'rV{r)e-'''-'' , (3.12a) 
^ (3.126) 



87r2A2 A;2 + A-2 ' 

(iii) J7d[ri,ri] is the disorder contribution to the total dynamic functional, in which 
we again transformed to Fourier k space, 

Jd[r^,n] = -DJ2j dtj^ dz j^{k-fi)U{k,z)e'''-^^ . (3.13) 

The Jacobian function J[fi\ in Eq. (|3.3|) depends on the discretization scheme of the 
path integral. It is simplest to adopt the causal discretization procedure in which the 
Jacobian is a constant, independent of and can be omitted, provided that simultaneously 
the ambiguous equal time correlator is defined to vanish, (see Ref.|29| for the details). 

{ft{z^)r]{z',t))=Q. (3.14) 

The averages of physical observables as well as the correlation functions of these 
observables can now be expressed in terms of functional derivatives of the generating 
function with respect to external fields that couple to these observables. For example, 
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n-point correlation function of an observable 0{fi, r^) can be obtained from the generating 
function Z[h{z,t)], 



Z[h{z,t)] = j Vf,{z,t) j Vr,{z,t) 



X exp 



-Jid[ri,ri] + I dt J dz O [fi{z,t),ri{z,t)) h{z,t) 



(3.15) 



by functionally differentiating n times with respect to h{z,t), 



O [fi{zi,ti),ri{zi,ti) ] ...O {riizn,tn),ri{zn,tr,) 

S''Z[h{z,t)] 



5h{zi, ti) . . . 5h{zm tj 



(3.16) 



{7i(zi,t,)}=0 

Hence, all the information about the dynamics of the flux line liquid is encoded in the 
generating functional which we will study in the hydrodynamic limit, below. 



4. Derivation of Hydrodynamic Description 

In the neutron scattering experiments, the neutrons interact with the magnetic fleld 
of the flux lines. The resulting scattering intensity is therefore proportional to the Fourier 



transform of the 2-point correlation function of the magnetic fleld. [31| |3C] The 

London equation relates the Fourier transform of the magnetic fleld components, along 
and perpendicular to z, to the flux line number and tangent density n{k,q,t), T{k,q,i), 
respectively, [|23| , p!8[1 



B.{k,q,t) = n{k,q,t) , (4.1a) 

1 + \\Mi{k^ + q2) 

BUk.q.t) = „^ P;[,{k)n{k,q,t) 

+ ^ ^t'n2 PaMnik.q.t) , (4.16) 

1 + X]_Mi{k^ + q2) 
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where -P^(A;) = Sab — kakb/k'^, P^bik) = kakb/k'^, and with the densities related to the flux 
hne position, 

N 

n{f,z,t) = '^6'^^^ {f-f,{z,t)) , (4.2a) 
1=1 

f(r,z,t) = 5^^5(2)(r-rl(z,t)) , (4.26) 

i=l 

Hence inelastic neutron scattering experiments are a direct probe of the dynamic structure 
function of interacting flux lines, i.e. of the flux line density-density correlation function. 

The transport coefficients like flux-line liquid viscosity and friction can also be ex- 
tracted from the dynamic structure function of these line liquids, and hence can be com- 



pared with the experiments that measure resistance. |T^ The interacting dynamic structure 
function also enters into the perturbative calculation of disorder corrected flux flow velocity 
in the high velocity regime. It has been argued that this relation holds beyond its expected 
regime of validity, down to low flow velocities. The dynamic structure function is therefore 
directly related to the current-voltage curves. []33|j3^ 

Because it is the flux-line density n(r, z,t) rather than the microscopic conformation 
fleld fi{z,t) that is measured in most of the experiments, it is more natural to work with 
the hydrodynamic density variables. By starting with the microscopic theory deflned by 
Jid[fi,'ri\i Eqs.(p?8| ) - (|3.13|) , and integrating out the microscopic conformational degrees of 
freedom ri{z, t) and fi{z, t), with the constraint that the densities n(r, z, t) and correspond- 
ing response fleld h{r,z,t) remain flxed, we derive the dynamic functional of the density 
flelds. With this new form of the dynamic theory the density correlation and response 
functions are easily computable and the approximations that are necessary can be more 
easily physically motivated, because the density and its correlations are directly observable. 

We begin by introducing a new density response fleld, n(r, z,t), 

N 

nik, z, t) = ksTD J2 [k ■ n{z, t) g-^^'^^^^^'*^ , (4.3) 

i=l 
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in addition to the physical density field n(r, z, t) already introduced in Eq.( [4.2a| ). The new 



response field n(r, t) will earn its name by generating dynamic density response Sanctions 
(see below and Appendix C). Here we will only study correlation and response functions 
of the number density field and therefore will trace over the tangent density field r(r, z, t) 
that is related to the fiuctuations of the magnetic field in the a6-plane, Eq. (|4.2^) . 

We reexpress the interaction parts of the dynamic functional, Eqs.( |3.1I| ), ( ^.131 ), in 
terms of the density fields, 

TV ^8 
jS,n]=iD / / J2^'^^Hr-niz,t))J2^j-—6'~^\r'-fjiz,t))V{r-r'), 

Jt.z Jr.r' „_i -^^ OTj 

(4.4a) 

n{f,z,t)fi{f',z,t)V{f—r'), (4.46) 



1 




t,z J r J r 



and the disorder contribution, 

f ^ Q 

Jd[r^,n]=iD Y,^^ 3^''\r-r,{z,t))U{r,z), (4.5a) 

Jt,z,r^^^ OTi 

= -r^ [ fi{r,z,t)U{r,z) . (4.56) 

f^B-L Jt,z,f 

It is important to remember that in the above equations h{r,z,t) and n{r,z,t) are not 
fields independent of the fundamental microscopic fields ri{z,t) and fi{z,t) and at this 
point are only a notational convenience. 

Although the mean-field theory of the flux-line liquid hydrodynamics is described by 
the average densities uq = {n{f,z,t)) and ho = {fi{f,z,t)), it is the fluctuations about 
the mean density that are nontrivial and are probed by the scattering experiments. We 
therefore introduce flelds, 

p{f,z,t) =n{f,z,t) -riQ , (4.6a) 
p(r, z, t) = h{f, z, t) — ho . (4.66) 

describing the deviations from the average uniform densities no = {n{r,z,t)) and no = 
(n(f, z,t)). 
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It is convenient to work in the canonical ensemble in which the number of flux lines 
is flxed (i.e. total magnetic field through the sample is constant), and for fixed sample area 
A, no = N/A. Substituting Eqs.( |4.6a7^ ) into Eqs.(4.4 ), (4.5 ), using the conservation of 
fiux lines identities, 

[p{f,z,t)=pik = 0,z,t)=0, (4.7a) 

J r 

p{f,z,t) = p{k = 0,z,t) =0 , (4.76) 



and ignoring the remaining constant terms that do not effect the dynamics, we obtain 
interaction part of dynamic functional in Fourier space, 

Ji[ri,ri] = I p{k,q,u)V{k)p{-k,-q,-u) , (4.8a) 
^ ^yi~^ I Pa{k,q,u)Vo,p{k)pp{-k,-q,-u) , (4.86) 

^Kb-L Jk.a.u; 



where we define, 



P= (P.P) , (4.9a) 

\ ^^'^) , (4.96) 

For the disorder contribution to the dynamic functional we obtain, 

Jd[n,ri] = tKt; L p{k,q,u)U{-k,-q)5{u) , (4.10a) 

Paik, q, u)ua{-k, -g, -u) , (4.106) 



1 

ksT Jk^q^oj 

where we defined disorder vector, 

u{k,q,u) = (U{k,q,u)5{u),0^ . (4.11) 

As discussed in Sec. 3, we construct the dynamic generating functional for the compu- 
tation of density correlation and response function, by choosing the operator Oa.{fi,ri) ~ 
Pa(f, z, i) and coupling an external field ha(f, z, t) to this density observable. 



Zd[h{r. z, t)] = fvR exp -J^d[R{z, t)] + j p{r, z, t) ■ h{r, z, t) 



(4.12) 
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In the dynamic functional, expressed in terms of densities, the disorder field, that couples 
linearly to the density field, can be integrated out exactly. It is important to note that in the 
dynamic formulation presented here we can integrate over the quenched disorder directly 
at the level of the dynamic generating functional Zd[h\. This is to be contrasted with 
static calculations where one must first compute the physical observables like free energy 
or correlation functions for particular realization of disorder and then perform the average 
over the disorder, (see Appendix A) This leads to the usual problems of disorder-averaging 
In or Z'^^ (since the static averages have to be normalized by the partition function Zd), 
problem usually handled using the "replica trick" . In the dynamic calculations presented 
here no such complications arise because the "dynamic partition function" Zd[h = 0] = 1, a, 
constraint enforced by the MSR Jacobian in Eq. (|3.3|) . We then integrate out the quenched 
disorder exactly, using the assumption of Gaussian disorder, Eq.(2.5 ). Upon averaging 
2d[h] over the quenched pinning potential -u, we obtain 



Z[h{f,z,t)]=Zd[h] , (4.13a) 



VR exp 



-J[R{z, t)] + / p(f, z, t) ■ h{ f, z, t) 



(4.136) 



where, 

J[R{z, t)]=Jo [R] + J.r^t \A , (4. 14) 



and 



with 



Jint[p\ = ttTTT; I pc{k,q,u)Kc,j3{k,q,u)pp{-k,-q,-u) , (4.15) 



''f') . (4.16) 

The interaction and disorder contributions have now been expressed as a quadratic 
function of the density fields pa- However, the computation of the dynamic generating 
function Z[h\ is still nontrivial because (i) the independent variables are R and not p^, and 
(ii) the noninteracting part of the dynamic functional Jq [R\ cannot be trivially rewritten in 
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terms of the density fields. To make progress however we can proceed with an uncontrolled 
variational approximation. We replace [i?] by an Ansatz that is Gaussian in the density 
fields 

Jo[-R] Jo[Pa] = 77—1 Pa{k, q,oj)T~Q{k, q,uj)pi3{-k, -q, -uj) , (4.17) 

^"■0 Jk,q,uj 

— * 

and use a measure Vp^ instead of VR in the functional integral. The matrix Fq,^ is 
determined by requiring that the correlation functions of the Gaussian Ansatz agree 
with the original noninteracting theory where the averages are performed with measure 
T>Rexp{—Jo[R]). The simplest requirement is that the structure functions in two theories 
agree. This completely fixes Fq,^ = S^^, where S^^ is the structure function for a single 
flux line. Combining with J'^^t this approximation then gives the full dynamic functional, 



^[Pa] = ni / Pa{k, q, uj) \noKai3{k, q, u) + ksTS^g ^{k,q,u;)\p/3{-k,-q,-uj) , 

ZKbI no Jk,q,u^ L J 

(4.18) 

from which all the dynamic correlation and response functions of density fields can be 
computed. 

In the next section we will treat the full dynamic functional more rigorously. We will 

show how the above ad hoc approximation emerges as a result of truncation of a systematic 
expansion of the dynamic functional in the density fields. 



5. Decoupling Flux Line Dynamics (Method of Auxiliary Random Fields) 

As we have already noted the density fields p(r, z, t) and p(r, t) are not independent 
of the fundamental microscopic fields fi{z,t) and fi{z,t) and at this point are only a 
notational convenience. The computation of the generating function in the above equation 
is made nontrivial by the highly nonlinear flux line interactions V{r), when expressed 
in terms of the kinetic flelds fi{z,t) and fi{z,t). To obtain interaction corrections to 
the noninteracting lines dynamics via a straight perturbation theory in powers of V{r) is 

19 



possible but is a nontrivial exercise. [3^| We note however that this interaction is simply 
quadratic when expressed in terms of the density fields. The idea is then to transform 
the functional integral from the microscopic variables to independent density fields. This 
can be accomplished using the method of auxiliary random fields that has been previously 
applied to treat both statics and dynamics of polymers and is described below. |l35[] 

We introduce a set of transformations that transform the functional of R{r, t) into 
a functional for the density fields. We accomplish this by introducing an independent 
auxiliary density fields ?/'(r, z, t) constrained by the functional 5-function to equal to the 
physical density fields p(r, z, t) for all r, z, t, 

Z[h{f,z,t)] = J VRe-^-^^^-^'-'^P^^^ j PV^JJ 5(iA(f,^,t) -p(f,;z,t)) (5.1a) 

r,z,t 

2)^g-^.n*[^1+/V^-^^^(^_^^ , (5.16) 

where the average is performed with the Gaussian dynamic measure e~^°[^] of noninter- 
acting fiux lines, with j7o[-R] given by Eq.( p.9| ). We observe that with this transformation 
the troublesome interaction is indeed quadratic in ip{r, z, t) and therefore in this sense 
can be treated nonperturbatively. The nontrivial part of the calculation reduces to the 
computation of the average of the functional (5- function with the Gaussian measure e~'^°^^\ 
Using the functional representation of the 5-function we obtain, 

(^d{^f-p)) = (^jv$e'I^-^'^-P'^^ (5.2a) 
= J V$e' I (e-' I ^-P^ . (5.26) 

Noting that p(r, z, t) = J2iLi Pii'^^ with pi{f, z, t) as the single line density, and since 
the average for each line is identical, Eq.(5.2 )further reduces to, 

(^Si^P-p)^ = J V$e'I'^-^[(^e-'I'^-^'') . (5.3) 
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Combining this result with interaction contribution to the hydrodynamic functional 
in Eq.(5.1 ) we obtain the dynamic generating function 



:[h] - j-Di^j 



J[ii;,cP]+ / ij{r,z,t)-h{f,z,t) 



(5.4) 



— * — * 

in terms of the hydrodynamic functional J'[ip., (f)], 

J[iJ, = J^ntW + n„r[^] - z y ^> , (5.5) 

with the new contribution 7ior[^] coming from the average of the 5-function in Eq. (|5.3| ) 
and is expressed in terms of a single flux line cumulant expansion, 

+ -^^ J <Pc.(pp(p^<PsT^^p^s + . . . , (5.6) 

where ra7a2 • is a connected correlation function of a single line density fluctuations 
and in coordinate space x = (r, z,t) is 

rK2...a™(^l'^2,---,a;m) = A{p^^{xi)pc,^{x2) . . . Pc,^{x^))o , (5.7) 

where the area factor A was inserted for convenience. In above we have also dropped the 
single line label 1 on pi since for a while we will be dealing with a single line density 
and no confusion should arise. The cumulant functions r^a2...am{xii^2, ■ ■■■,Xm) can be 
systematically computed for any m (see Appendices B and C). 

We note that the density field -0 appears only linearly and quadratically in Eq. (|5.4| ) 



and therefore can be integrated out exactly. This leads to the dynamic functional that is 
expressed not in terms of the physical densities tj) but in terms of the auxiliary fields (f), 
conjugate to ip. Since the hydrodynamic functional is derived in the presence of sources h 
the density correlation functions and response functions of can still be easily computed 
as functional derivatives with respect to these sources (see Eq. (|3.16|) ). Upon integrating 
out if: we obtain, 

Z[h] = j V^e-^^^^'^y (5.8) 
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where, 

J{4>, h\ = - y i'<Pc. + h^)K-^{i<Pp + hp) . (5.9) 

The above equations define the hydrodynamic field theory of fiux fines derived directly 
from the microscopic interacting single line dynamics. No approximations have been made 
up to now, and both the fiux line interactions and interaction with disorder have been 
treated exactly. To fully solve the hydrodynamic theory we need to compute the generating 
function Z[h] by integrating over (f). Since the resulting theory has a form of an interacting 
field theory, it cannot be solved exactly. However, one could try to treat the nonlinear 
interactions perturbatively with various techniques. In the long wavelength limit (which 
can in principle be studied with sufficiently low angle neutron scattering experiments) the 
(f) fiuctuations can be treated with renormalization group methods as was done for the 



statics of polymers solutions by Ohta and Nakanishi||35[] . This procedure would lead to 
renormalized vertex functions, T^'^\ corrected by the thermal fiuctuations of (f). As will be 
shown in the next section, some of the fiux line fiuctuations and interactions are nontrivially 
taken into account even if we truncate the hydrodynamic functional at the quadratic order. 
The resulting truncated theory can be solved exactly, allowing calculations of dynamic and 
static properties of the interacting fiux-line liquid. 



6. Gaussian Approximation to the Dynamic Functional 

In the previous section we derived the expression for the dynamic generating func- 
tional. It is expressed in terms of an infinite series of interactions in the auxiliary field 
cj) and in principle allows for calculation of any hydrodynamic correlation function of the 
fiux-line liquid. Near Hd, where the magnetic field is weak and the fiux-line liquid is 
dilute, the fiuctuations in the density fields are large, i.e. comparable to the average line 
density. In this critical region the nonlinear interactions must be carefully taken into ac- 
count. Using renormalization group methods the dynamics of dilute line-liquids can be 
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studied, as was done for the statics. Similarly, near Hc2i where the flux-line liquid is 
dense, the vortex cores begin to overlap, the effects of nonlinear interactions are large and 
renormalization group treatment is again needed to take the nonlinearities into account. 
13^ Here, however, we study the dynamics of a semidilute flux-line liquid, away from such 
critical regions, i.e. in the regime where Hd « H « Hc2- In this regime uq is large, 
hence the method of steepest descents applied to Eqs.( ^78D ,( |5.9|) allows us to treat Z[h] in 
mean-field theory. The nonlinear interactions can then be treated in perturbation theory, 
but in this semidilute regime we do not expect qualitative corrections to the results of our 
mean-field approximation. In this section we will therefore use a Gaussian approximation 
for computation of Z[h] by truncating expansion of the dynamic functional Eq. (|5.6|) , at 
quadratic order in (p. 

The remaining T^p{k,q,ui) vertex function has a clear physical interpretation of a 
matrix whose (2,2) and (2,1) components are the noninter acting dynamic structure func- 
tion S°{k,q,uj) and response function S°{k,q,uj) of the fiux line liquid, as can be seen 
from Eq.(|5.7|) . V)^^{k,q,uj) is computed and analyzed in various regimes in Appendix C. 
With the truncation at quadratic order we obtain the interacting stucture and response 

~* ~* (2) ~* 

functions Sa/3{k, q,uj) in terms of the noninteracting ones S^i^{k,q,uj) = r^^(/c, a;) and 
therefore will be able to study the effects of fiux-line interactions and quenched disorder 
on the dynamics of the fiux-line liquid. 

In the Gaussian approximation it is more convenient to return to the Eqn.f ^T^) and 
integrate out the cj) field, thereby producing a hydrodynamic functional jT'fV'] expressed 
directly in terms of physical density fields ?/'. Upon performing all the calculations in 
Fourier space {q and k are reserved for wavevectors in the z and the transverse directions, 
respectively) we obtain. 



Z[h] = J Dtf expjy^ {^]^i;oc{k,q,uj)k^p{k,q,uj)'ipp{-k,-q,-u) + 

+'i(;o,ik,q,uj)h^{-k,-q,-uj) } . (6.1) 
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where, 



K^p{k,q,u) = T^Kc,p{k,q,u) + —{S°)J^{k,q,u^) . (6.2) 



An integration over ij) leads to the dynamic generating function within the Gaussian ap- 
proximation, 

Z[h]=e^Y>\\ I ho,{k,q,uj)k~p{k,q,u)hp{-k,~q,-uj)\ . (6.3) 

L ^ Jk,q,w ) 

FunctionaUy differentiating twice with respect to external field h{k,q,uj), we obtain the 
interacting correlation/response function matrix Sap{k, q, uj) = K~p(k, q, w), 

\b[-k,-q,~uj) b[k,q,uj)J 

where ^22 = S{k,q,uj) and 5'i2 = S{k,q,uj) are the interacting hydrodynamic structure 
and response functions, respectively, for the flux line liquid in the presence of disorder, 

^.r . F{k, q)S{u;)\S"{k, q, u^)\''nl/ [ksTf + UpS^k, q, uj) 

b{k,q,uj) = ^—3 — , (6.5a) 

\Vik)S"{k,q,u)nJkBT + l\^ 

^it , V{k)nl\~S"% g, c^) | Vfc^T + no~S"{k, g, u) 

b(k,q,uj) = =r-; — =; . (6.00) 

\V{k)S"ik,q,Lo)no/kBT + l\^ 



Note that in Eq. (|6.4|) the Sn component vanishes, {tl^lf, z,t')tlj{f ' , z' ,t')) = 0. This 
is a consequence of the fluctuation dissipation theorem and causality, and implies that 
dynamics of the density field 'i/'2 encoded in the effective dynamic functional in Eq. (|6.3| ) 
can equivalently be described by a linear differential equation for the density field "0. 



7. Details of the Dynamic Structure Function 

In this section we will analyze the dynamic structure function S{k, g, w). Information 
about the behavior of S{k,q,uj) can be obtained from S{k,q,uj) by using Fluctuation 
Dissipation Theorem (FDT), 

dtV^S ik, z, t) = rg (k, z, t) - rg (fc, z, -t) . (7.1) 
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The interacting dynamic structure function in Eq. (|6.5a|) is expressed in terms of the non- 
interacting structure and response functions, S°{k,q,uj),S°{k,q,uj) of the vortex Hquid 
and depends on the correlation function of the disorder F{k, q) as weU as the interhne 
interaction V{k). 

The details of a single line dynamics are presented in Appendix C. There are two 
regimes with very different dynamic behavior that are seperated by a characteristic time, 
i Rouse = L'^/{De). The Rouse time is the time required for the center of mass of the 
flux line to diffuse its transverse radius of gyration. Equivalently, for times larger than 
tRouse the diffusion of the center of mass dominates over the internal mode dynamics. The 
dynamics of a single flux line for t ^ tRouse and t <ti tRouse is depicted in flg. ^. 

We find that for times larger than the Rouse time, tRouse = L'^/{De) the diffusion 
is dominated by the center of mass mode. In this case the average transverse distance 
diffused is simply 



r^W = f^{[ro{t)-fom') = f-^^^t'/- . (7.2) 

We note that the flux line diffusion coefficient DksT/L is scaled down by L with respect 
to the point vortex diffusion constant OksT in the a6-plane. This 1/L behavior has been 



previously derived in the appendix of Ref. in terms of a simple model of point vortices 
coupled in the z-direction. 

In the opposite limit, t < tRouse, the flux line diffusion is dominated by the internal 
modes. For two points on the flux line separated by a distance z <^ (Det)^/^ the auto- 
correlation function becomes indendent of z. This allows us again to extract a transverse 
diffusion distance 



rhit) = ]ll{[r{z.t)-f {0,0)]') ^ ^2kBT (^t^ cxt'/\ (7.3) 

but now with an anomalous t-dependence. This behavior can be understood as two "ran- 
dom walks" on top of each other; the flux line segment executes a random walk on the flux 
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line conformation, which can be thought of being generated by a random walk (fictitious 
dynamics along the 2;-axis). For z ^ (Det)^^'^ one gets "diffusion" in the time-like variable 



z. The corresponding square root of the mean square displacement y ^^-^^-^1^1 corresponds 
to the projected 2D radius of gyration of a flux line of length 1^1 . At intermediate scales 
z > (Dety/'^ the 2;, t-dependence of the segment-correlation function has to be taken into 
account, and one flnds 



^2, 4A;bT, , , fDe\t\ 



([r(z,t) -r(0,0)r) = ^kl/ ( ^ ) , (7.4) 



with f{x) shown in flg. ^ 

The dynamic crossover described above can be equivalently extracted from the be- 
havior of the LV-pole in the noninteracting structure function. For k^L <^ 1 the center 
of mass diffusion dominates the hydrodynamics and the noninteracting relaxation rate is 
^cm{k) = DhsTk^ / L. In the opposite limit of k^L ^ 1 the internal modes dominate 
the dynamics and lead to the relaxation rate F^^ = {kBT)'^Dk^ / (4e). This internal-mode re- 
laxational rate can be equivalently rewritten in the form of the center-of-mass relaxational 
rate F^ = DksTk^ /L{k), with an effective fliix line length L{k) = j^^- Not surpris- 
ingly, L{k) is approximately the length in the 2;-direction corresponding to the xy-length 
scale of k~^. The crossover from a dynamics dominated by the center of mass motion to 
a dynamics governed by the internal modes (or equivalently from L to L{k)) occurs when 
the wavelength (transverse length scale) ~ becomes smaller than the static transverse 
"diffusion" length ^j2kBTL / e (the 2D projected radius of gyration). 

Substituting the expressions for the noninteracting correlation and response functions 
in these two regimes from Appendix C into Eq. (|6.5a| ), we obtain the dynamic interacting 



structure function in the Gaussian approximation. 
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7.1. Center of Mass Dominated Regime 



In the wavevector regime, k L ^ 1 , where the center of mass mode dominates 
we obtain 



b[k,q = [J,u))^ ^ -\ 2 • ('-5) 

u;^ + [DksTk- {i + ^)] + 

Concentrating on the first term which describes the dynamics in the absence of disorder we 
observe that the center of mass, noninteracting relaxation rate rcm(^) has been additively 
renormahzed by the interaction between the flux lines. We define an interaction length 
scale in the 2— direction, Lj{k) = {uoV (k) / kBT)~^ , and note that the crossover between 
the noninteracting and interacting dynamics occurs when L > Li{k). (see fig. |]) 

Physically, Ljik) is the fiux line length (or equivalently sample thickness) beyond 
which there are multiple interactions between the flux lines. For samples thicker than 
Li{k) (if the crossing barriers are large) signiflcant line entanglement will take place. This 
length is therefore analogous to the entanglement length discussed in Ref.|l| . 

The expression for L/(/c) derived in the Gaussian approximation will be corrected 
by the higher order interactions appearing in the dynamic functional expansion, Eq. (|5.6|) . 
The flux line dilute regime can be treated using renormalization group as was flrst done for 
the statics in Refs.p],|2B|. It is found that in the dilute limit the interacting theory becomes 
asymptotically free, with the flux line interactions renormalizing to zero, logarithmically 
with length scale. Following the renormalization group equations down into the dense 
regime and matching, gives the effective interaction V ~ 47r(/c5T)^/e/ ln(l/7ioA^), which 
is independent of the original bare interaction. In the dilute regime (dropping the unimpor- 
tant logarithmic factor and constants of order unity) we therefore obtain L/ ^ e/{2kBTno). 
This length is exactly the entanglement length deflned by the static transverse wandering 
being on the order of the average line spacing. 
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For three dimensional samples (L oo), the flux line density relaxations are dom- 
inated by the interactions with the rela:xation rate given by the rate for flux lines with 
length Li{k), 

rf„ ^ ^ . (7.6a) 

Li{k) 

~ noV{k)Dk'^ , in the dense regime , (7.66) 

~ , m the dilute regime . [J .be) 

We thus find that while for the noninteracting fiux line the center of mass diffuses 
very slowly, with the rate vanishing as l/L, the relaxations in the interacting line liquid 
are independent of the flux line length L. In the interacting theory the 1/L dependence 
of the diffusion coefficient gets cutoff by the interaction (entanglement) length L/. Within 
the Gaussian approximation we therefore find that the fiux-line interactions speed up the 
dynamics of the line liquid, probably because they lead to a stiffer response to density 
inhomogeneities . 

By comparing with the results of the static structure function obtained with the 
Gaussian approximation (see Appendix A and the second part of this section) we find that 
fiux-line liquid kinetic coefficients are not modified and the renormalization of the statics 
is solely responsible for the increase in the relaxational rate Fcm- Therefore, within the 
Gaussian approximation no viscosity is generated. We expect, however, that this will be 
corrected by higher order interactions, which also become important in the dense limit, 
near Hc^- It should be possible to control these interactions via renormalization group 
or mode coupling theory. If the line crossing barriers remain high, the additional effects 
of fiux line entanglement, not taken into account by the Gaussian approximation, will 
play a major role in the dynamics and probably can be described by these higher order 
nonlinear ities. We expect that entanglement and caging effects to become important in 
the regime of high line densities and large fiux line crossing barriers, which is possible 
at intermediate field strengths. These effects should have a slowing down effect on the 
dynamics, and will compete and eventually swamp the relaxational rate increase found in 
Eq.(^) . 
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1.2. Internal Modes Dominated Regime 

In the regime, k^L^ 1, [37| where the internal modes dominate we find 

2noCx{k, q) 



S{k,q,uj) 



+ 



uj^P^{k,q)+^^{Kq) 



keT 



5{u)Fik,qy 



A{k)K{k,q) 



1 + 



A{k)K{k,q) 



2 ' 



(7.7) 



where we have introduced 



a{k, q) = A{k)^^Vr 



k ' 



P^k,q) 




1 + ^^(*)6„ 

ksT 



bs-fHk.q) binoV{k)A{k) 



(7.8a) 
(7.86) 



+ 



h ri bs 

noV{k)A{k) 



keT 



ksT 
(26062 - bl) 



262+ 



(7.8c) 



and, A{k) = Ae/{kBT^k^) and = DiksTyk'^ /{4e). The coefficients 6^ are calculated 
in Appendix C and for 2eq/ {ksTk'^) are constants. The interacting dynamic structure 
function above consists of two parts. The first term is the thermal contribution to the 
density-density correlation function and has the standard Lorenzian shape. We note that 
the coefficient in front of the w^-term depends nonanalytically on wavevectors k and q 
and is therefore a breakdown of traditional hydrodynamics, here due to the effects of the 
internal modes. 

The second term in Eq.( |7.7| ) is the elastic contribution to the dynamic structure func- 
tion arising from quenched pinning disorder. This contribution leads to a time independent 
persistent contribution to the structure function. While the density correlations of the un- 
pinned fraction of lines decay in time as the lines move around, the pinned density fraction 
has correlations that are time independent and have spatial correlations of the random 
potential. 
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The linewidth 7(/c) has the exphcit form 



7(^) = (^noV{k)De + ^-^Dk' ] (7.9) 



Note that the /c'^-term is independent of the interaction between the flux hnes and is there- 
fore just related to the internal dynamics of a single flux line. We find that no interaction- 
generated viscosity appears. However, flux-line interaction modifies the relaxation rate of 
a single line in an important way. 'j{k) describes the crossover of the relaxation rate from 
the [kBT^Dk^ / (4e) behavior of a noninteracting line liquid to noV{k)Dk'^ relaxation rate 
of an interacting liquid. This crossover is similar to the Bogoliubov crossover in statics [|l|]. 
The interacting dynamics is again diffusive and aside from some constants is similar to the 
center of mass dominated regime considered in beginning of this section. 

7.3. Static Structure Function 

The static structure function Ss{k,q) is an equal-time density correlation function, 
and therefore can be obtained from the dynamic structure function, 

Ss{k,q) = S{k,q,t = 0) , (7.10a) 
^S{k,q,u). (7.106) 

Applying these equations we obtain the interacting static structure function of the vortex 
liquid in the presence of disorder (see also Appendix A), 

S.Ck. - + ^ ( V , (7.11a) 

l + noV{k)S^{k,q)/kBT {kBT)^ \l + noV{k)S",{k,q)/kBT J 

+ Fik, g) ( ^ ] ' , (7.116) 



g2 + {qB{k)/kBTY ' + {qB{k)/kBT) 

where qsik) is the Bogoliubov spectrum of the corresponding bosons, |jl| 



(lB{k) 
kBT 



kBTk'^\'^ ^ noV{k)k^ 



2e 



1/2 



(7.12) 



30 



The equations for the structure functions that we derived here are vahd for general 
type of disorder, characterized by disorder correlation function F{k,q). The physically 
relevant cases are: 

(i) Point disorder due to oxygen vacancies (uncorrelated) 

F(f, z) = AoS'^^\r)6{z) , (7.13a) 
F{k,q) = Ao. (7.136) 

(ii) Columnar defects (line-correlated along z-axis) 

F{f,z) = AiS'^^^r) , (7.14a) 
F{k,q) = AiS{q) . (7.146) 

(iii) Grain/twin boundaries (plane-correlated, with normal n in xy-plane) 

F(f,z) = A25(f • n) , (7.15a) 
F{k,q) = A2S{z-{kxn))6{q) . (7.156) 

Eq. (|7.7| ) also depends on the interline interaction V{k) and is valid for a general range 
of interactions or equivalently for arbitrary line density. In the limits of low and high line 
densities V{k) reduces, respectively, to 

Vik)^^, A;A<1, (7.16a) 

n^)^^;^ = ^o, kX»l. (7.166) 

Eq.( [7.16a| ) is valid for H ^ Hd where the lines are much farther apart than A, and 
Eq.( [7.16"^ holds for H » Hd where the average interline distance is smaller than A. 

Specializing our general result for the interacting static structure function, Eq.(7.11 ) 
to point disorder and to short-range line interaction, we recover the result of Nelson and 



LeDoussal, 23 



Sl^{U)= ""^^^^/^ -fAo ""-^ I , (7.17) 

q^ + {qB{k)/kBTy \q^ + {qBik)/kBT)\' 

They utilized the boson mapping to obtain this static result. For pure superconductors 
this result reduces to the original result of Nelson and Seung|^] also obtained via the boson 
mapping. 
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8. Phenomenological Hydrodynamics 

As explained in the introduction, one can also take a more macroscopic, phenomeno- 



logical approach to derive equations of motion for the hydrodynamic density fields. [18| 
We take this approach in this section with the intent to subsequently compare the results 
of our kinetic theory derived in previous sections with the hydrodynamic approach. 

As already emphasized in Sec.^ , on long time and space scales the important, slow de- 
grees of freedom that characterize the flux-line liquid are the number and tangent densities 
of the vortex liquid, 

N 

n{f,z,t) ='^6'^^^ {f-f,{z,t)) , (8.1a) 

i=l 
N 

t{r, z,t) = J2 dzn{z, (f - n{z, t)) , (8.16) 
1=1 

together with the current fields (r, z, t) = UoVaif, z, t), jai,{r, z, t) 

N 

j:ir,z,t) = Y,dtna{z,t)6^^Hr-n{z,t)) , (8.2a) 

i=l 
N 

j*fe(f, z,t) = J2 [^tr^a{z. t)d,r,b{z, t) - dtnb{z, t)d,na{z, t)] 5(2) (f - n{z, t)) (8.26) 
1=1 

which transport the number and tangent densities n(r, z, t), taif, z, t), respectively. As for 
the hydrodynamics of a liquid of point particles, the dynamics of a line-liquid is constrained 
by dynamic continuity equations arising from local conservation of the density fields, 

dtn{r, z, t) + aj^(f, 2, t) = , (8.3a) 
dtta[f, z, t) + d^ii^ir, t) - dj^ir, z,t)=0. (8.36) 

Furthermore, n{f,z,t) and ta{r,z,t) are not completely independent. Continuity of flux 
lines introduces a spatial constraint between the number density and the longitudinal part 
of the tangent density, 

dMr,z,t) + data{r,z,t) = , (8.4) 
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which in the language of bosons plays the role of a temporal continuity equation for the 
conservation of the boson density along imaginary "time" z. Eqs.(8.3 ), ( p.4| ) can be easily 
verified by substitution using the microscopic definitions in Eqs.(8.1 ), (8.2.) The dynamics 
of a fiux-line liquid is therefore governed by only two (as opposed to three) hydro dynamic 
variables, the number density n(r, z, t) and the transverse part of the tangent density 
F(f,^,t). 

It is also enlightening to interpret the above continuity equations in terms of the un- 
derlying electromagnetic fields. At length scales much larger than the London penetration 
lengths, » X±, » A, Eqs.(4.1 ) (or equivalently, the long wavelength limit 
of London equation) lead to a simple relation between the magnetic field and the line 
densities, 

B{f, z, t) = z(f)on{f, z, t) + (f)ot{f, z, t) . (8.5) 

It is easy to see that Eqs.(8.3 ) are just the z and x-y components of Maxwell equation 
dtB/c + V X i? = 0, with the electric field E, which results from motion of fiux lines, 
related to the current fields of Eq.(8.2 ).[^ The spatial continuity equation, Eq.( ^.4|) , is 
equivalent to V ■ S = 0. 

A closed set of hydrodynamic equations is obtained when we supplement the continuity 
equations with the constitutive equations for the currents. The constitutive equations are a 
statement of Newton's 2nd law, cast in terms of the hydrodynamic variables. The relations 
equate the rate of change of the velocity to the forces acting on the fiux-line liquid. These 
forces include the frictional forces, due to fiux lines' interaction with the underlying lattice 
and weak disorder, and the pressure gradients due to nonuniformity in the density of the 
vortex liquid. For simplicity we will treat the case of small Hall angle and therefore neglect 
the component of the velocity response perpendicular to the forces. At long times, and for 
large frictional forces the velocity fields quickly decay to their steady-state value and the 
inertia term can be ignored. The resulting time-independent constitutive relation equates 
the frictional forces, proportional to velocity fields, to the pressure forces, expressed in 
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terms of the density fields. A simplified equation for the velocity of the number density 
current is given by, 

h~V^l-Vzd^^)Va{r,Z,t) = -Uoda^— T^+riodzj—; - + f;^^\f,Z,t)+Ua{f,Z,t) . 

dn{r,z,t) dta{r,z,t) 

(8.6) 

The terms on the left hand side of the above equation represent the frictional and viscous 
forces acting on the flux lines, which are balanced by the pressure gradient forces, the 
external forces fext and a random noise force appearing on the right hand side of the 
equation, f^xt might include the Lorentz force, fext = —no4>oZ x due to the charge 
current je coupling to the magnetic field of the flux-line liquid. Here we are interested in 
the equilibrium regime, fext = 0. We take iya{r,z,t) to be a Gaussian zero-mean noise 
with the correlations determined by the Fluctuation Dissipation Theorem, 

{ua{k,q,uj)ub{k\q',uj')) = (8.7a) 
= ksT (7 + ryP + r],q^^ {2n)^S^^^ {k + k')5{q + q')5{uj + uj')5ab , 

{Pa(k,q,uj)) = Q. iS.lh) 

For the hydrodynamic description it is sufficient to assume that the effective Hamil- 
tonian, 7Y[n, t] is an expansion in powers of hydrodynamic variables n(r, z, t) and t(r, 2, t), 
which for simplicity we truncate at the quadratic order. It is more convenient to expand in 
terms of the fluctuations p(f, z) = n{f, z) —Uo and r(r, z) = t{f, z) — to around the average 
values Tio, to, so that the linear terms can be eliminated. Ignoring disorder (because it 
effects can be easily included) and assuming translational invariance, we obtain in Fourier 
space, 

n[p(k,q),f{k,q)]= [ l\K,{k,q)p^ + K2{Kq)r''] , (8.8) 

Jk,q ^ L J 

where functions Ki{k, q) and K2{k, q) are related to the static density correlation functions. 
Imposing the constraint of Eq.( p.4| ) in Fourier space, qp{k, q) = —k ■ r{k, q) to reexpress H 
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in terms of independent hydrodynamic variables, we easily compute the static structure 
function, {p{k, q)p{—k, —q)), and the tangent correlation function {Ta{k, q)Tb{—k, —q)). 

ksTk"^ 



S%Kq) 



k''K^{k,q) + q^K2{k, q) 



(8.9a) 
(8.96) 



This allows us to reexpress Eg. ( p.8| ) in terms of the static structure function S^{k,q) and 
the transverse part of the tangent density correlation function TT{k, q), 



K,{k. q) = Pll,{k)^^ + P^Mk)^S%k, q) 



'H[p{k,q),T{k,q)\ 



ksT 



k,q 



(^S:\k,q) - ^T-\k,q)^ p^ + T-\k,q)r' 



.10) 



Combining the above equation with Eqs.(|8.6|) , ( ^.4|) , (|8.3a|) we find that the dy- 
namic equations for n{k, q, t) and t{k, q, t) decouple. Here we will only concentrate on the 
hydrodynamics of n(/c, q^ t) which is governed by a simple Langevin equation. 



dtn + [nikBTD{k,q)eS:\k,q)]n = C{Kq,t) , 



l.ll] 



where we defined the hydrodynamic diffusion coefficient D~^{k,q) = 7 + r]k'^ + rj^q^ and 
the noise C{k,q,t) = inoD{k,q)k ■ V{k,q), which is also Gaussian with zero-mean and the 
correlations determined by those of V, Eq.(8.7 ), 



(C(fc, uo)Q{k\ q', uj')) = ksTDik, q)^ nl{2ix)H^''^ {k + k')5{q + q')5{u + uo') 
{ak,q,u;)) = 0. 



(8.12a) 
(8.126) 



These equations then lead to the hydrodynamic structure and response functions previously 



derived in Ref.l8 



2nlkBTk'^D{k,q) 



u;^ + (nlkBTkW{k,q)/Ss{k,q) 

nlkBTk^D{k,q) 
iu + n'lkBTk'^Dik, q)/ Ss(k, q) 

35 



2 ' 



(8.13a) 
(8.136) 



We are now in a position to compare the results of the phenomenological model of 
hydrodynamics with our kinetic theory of the flux-line liquid. Making the comparison in 
the physically most relevant regime dominated by center-of-mass motion, Eq.( [7.5|) , we 
find the following identifications, 

D < — > noD{k, q = 0) , (8.14a) 

[l ksT ) S,(k,q = 0) 

As we already observed in Sec.|^ no viscosity is generated within the Gaussian approxi- 
mation employed here and we are only able to establish a simple relations between the 
parameters of our kinetic model and the phenomenological theory, 

7^-^, (8.15a) 

Eq.( p.l5a|) is physically appealing in its identification of the phenomenological friction 
coefficient 7 with the inverse of the kinetic diffusion coefficient D. Eq. (|8.15^ ) leads to the 
the static structure function, in agreement with the result obtained from boson analogy. 
The above equations show that at least within the Gaussian approximation, in the center- 
of-mass dominated regime the dynamics is modified only through the statics, i.e. the 
kinetic coefficients are not renormalized in this order of approximation. 



9. Conclusions 

In this paper we have formulated the dynamic theory of the fiux line liquid phase 
directly from the kinetic theory of individual, interacting fiux lines. We have used the 
resulting theory to study the dynamic structure and response functions of the line liquid 
in the presence of various types of pinning disorder. In order to solve the theory we 
employed a Gaussian approximation to the dynamic functional, which should be valid 
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at intermediate flux line densities or fields H^^ <^ H <^ Hc2, where the effects of large 
fluctuations are not as important. 

We expressed the interacting dynamic structure function in terms of the single line 
structure functions. In the long time limit, t ^ tuouse and/or for transverse wavelengths 
larger than the line wandering, the center of mass mode dominates over the internal mode, 
and we recover the hydrodynamics of rigid rods. While for noninteracting lines the diffusion 
coefficient vanishes as L —>■ oo, the interactions between the lines lead to an increase in the 
relaxation rate of the line liquid. The diffusion rate remains flnite for any L, with L cutoff 
by the interaction length Li{k). In the dense limit Lj is determined by the interactions, 
with the two-body interaction giving only a first order estimate to this length. A detailed 
calculation that takes into account higher order interactions is required in this dense limit. 
In the dilute regime, the renormalization group calculations lead to a line interaction that 
vanishes as an inverse of a logarithm of the length scale and the microscopic interaction 
drops out. In this case L/ becomes just the entanglement length defined in Ref.|l|. In either 
case for, t ^ tuouse, our results are therefore in agreement with the phenomenological 
model of Marchetti and Nelson, but unfortunately the Gaussian approximation is not 
accurate enough to generate the viscosity of their model from our kinetic theory. 

In the opposite limit of short times and/or ^sTk l \ g^d that the internal 
modes dominate. In the presence of interactions we find that the /c^ relaxation rate of 
a noninteracting fiexible line liquid crosses over to a diffusive /c^ relaxation rate. In con- 
trast with the phenomenological model, however, a wavevector-dependent coefficient of 
the term is generated. This coefficient and by this constitutes a break down of 

conventional hydrodynamics due to the internal mode fiuctuations. 

We find that the Gaussian approximation successfully reproduces the static Bogoli- 
ubov "spectrum" previously obtained via boson analogy. For the statics this approximation 
is therefore equivalent to an infinite summation of all ladder graphs in the single line pic- 
ture, and then treating the resulting field theory in the mean-field approximation. The 
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dynamic results of this approximation have the structure of a dynamic Random Phase 
Approximation (RPA) . In addition to the effects discussed above this approximation leads 
to a disorder generated perturbative correction to the dynamic structure function, which 
agrees with the phenomenological approach. 

It is well known that preserving the discreteness of the flux lines is crucial for the 
correct treatment of the low temperature ordered phases like the vortex glass and Bose 
glass. It is likely that in the liquid phase, however, much of the dynamics can be described 
in terms of the density fields which coarse-grain over the discrete line coordinates as we 
have done here. In this description slow dynamics for high line crossing barriers and 
entanglement can in principle be incorporated by higher order nonlocal interactions in the 
densities, although only in some average sense. Also line crossing and recombination are 
present in the theory and can be controlled in an average sense by the strength of the 
repulsive flux-line interaction. This description of a flux-line liquid can be improved by 
also introducing a local tangent fleld that keeps track of the local orientation of the fliix 
lines as we do for the statics in Appendix A. 

We expect that the results derived in this paper will be corrected by the higher 
order interaction which can be taken into account using renormalization group or mode 
coupling methods. These corrections will undoubtly lead to renormalization of the kinetic 
coefficients. We also expect that the slow dynamics resulting from the lines being caged by 
its neighbors and from fliix line entanglement effects, can be reproduced by these higher 
order interactions. It is quite likely, however, that these effects will turn out to be non- 
perturbative. It is also possible that the k'^ correction to the diffusion coefficient will be 
absent even beyond the Gaussian approximation. The interactions might generate only 
nonanalytic terms in k and therefore the flux-line liquid viscosity will be strictly zero, 
leading to a breakdown of the phenomenological model described in Sec. 8. In any case we 
believe that the formulation of hydrodynamics in terms of kinetic theory of lines introduced 
in this paper will be useful for a more detailed understanding of flux-line liquid phase at 
and possibly away from equilibrium. 
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Appendix A. Calculation of the Interacting Static Structure Function 

In this Appendix we derive the interacting static structure function of Eq.(7.11 ) 
within the static equilibrium formulation, directly from the microscopic Hamiltonian 
Eq.( p.l| ), using methods similar to the ones used in the main text to obtain the dynamic 
structure function. 

We begin with the microscopic Hamiltonian, 



n{z),z) . (A.l) 



^ = ^E/ dzid.r.f + l [ dzV{n{z)-r,{z)) + J2 f dzU{ 

^ z=l ^ ^^J = l ^=l 

Analogously to the dynamics calculation we want to derive a macroscopic Hamiltonian in 
terms of number and tangent densities, p(r, z) = n{r, z) — uq and r(r, z) = t{f, z) — to, 
where, 



N 

n[r, z) = 

i=l 
N 



J2S^'Hr-r^{z)) , (A.2a) 



r(f, z) = Y, ^^^'^ (r- r,{z)) , (A.26) 

i=l 

The tangent density field can be written as a sum of its longitudinal and transverse parts, 
r(r, z) = T^{f, z) + r'^if, z) with r^(r, z) = -P^T&(r, z) and rj(r, z) = P^b'^bifi z), where 
^abi^) ~ ^ab — kakb/k? aud P^ij{k) = kakb/k"^, in Fourier space. 

As already mentioned in the derivation of dynamics, the longitudinal part of the 
tangent density field T^{f,z) is related to the number density n{r,z) by the "continuity 
equation" , 

—p{f,z)+Vr-f{f,z) = 0, (A.3) 

which is a statement that the flux lines do not end within the sample, or equivalently in 
terms of the magnetic field is Maxwell equation, V ■ -B(r, z) = 0. This equation is clearly 
trivial for r'^if, z), but gives a relation between p(r, z) and r^(r, z). Clearly then it is not 
necessary to keep track of the longitudinal part of the tangent density since its correlations 

40 



can be easily obtained from those of the number density using Eq. (|A.3| ) . We derive an 
effective Hamiltonian in terms of static density fields p(r, z) and r^if, z) by starting with 
the microscopic Hamiltonian Eq.( [A.l|) and integrating out microscopic degrees of freedom 
fi{z) as we did for the dynamic calculation. 

The microscopic Hamiltonian can be rewritten (aside from irrelevant constants) in the 
suggestive form in terms of the relevant density fields, 

where Ti.o[fi] is the Hamiltonian of noninteracting lines, the first term in Eq.( [A.l|) and. 



T-Cint[p,r]= p{f,z)p{r",z)V{f-r")+ p{f,z)U{f,z). (A. 5) 

J r,r ' ,z J r,z 

Since we are looking for a description in terms of the density fields, we construct a gen- 
erating functional by introducing external fields /i(r, z) and hrif, z) that couple to p(r, z) 
and T^{f,z), respectively, 

Zd[K K] = J Vfiiz) exp + j {p(^^ z)h{f, z) + f(f, z) ■ K{r, z)^ 

(A.6) 

We introduce unity inside above expression in the form of functional 5-functions, con- 
straining the auxiliary fields -(/^(r, z) and i^rif, z) to equal the physical densities fiuctuations, 
respectively, 

Zd[h, K] = j Vn{z)e-^'^"^^-P'^^-^'^P^^e^r^M''''^'^^^^^^ 

X J V^VtprYl5{tp{f,z) - p{f, z)) 5 (^iJr{r, z) -T^{r,z)^ 

X (^5{i;{f,z) - p{f,z))6 (^ifr{r,z) -T'^{f,z)'^^^ , (A.7) 
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The interactions and disorder are quadratic and linear functions of auxiliary density 
fields and the nontrivial part of the calculation reduces to averaging the functional 5- 
functions with the Boltzmann weight exp[— /37io] of noninteracting lines. These averages 
decouple to single line averages, 



6 {tpif, z) - p{f, z)) S{iprir, z) - r'^(f, z)) 

o 

-, N 



V(j)V(j)re 



i J ((j}ip+$T-^T) / -i J{<l>pi+(j}T -Tti 



(A.8) 



Inserting above equation into Eq.( [A.7| ) we obtain the partition function (for a fixed 
realization of disorder) in terms of an effective Hamiltonian expressed as an expansion in 
macroscopic density fields, 

Zd[h, hr] = J V^PVifr J V(j)V(fr CXp -Hf?/', ifr, </>, ^r] + J {tph + tfr ' K) , (A.9) 

where the effective Hamiltonian is given by, 

n[ij, i^r, (t>, $r] = n^ntii^] + noTlc/}, $r]+i j (</'^ + " ^r) , (A.IO) 

with the new contribution noT[(f), (f)^] coming from the average of the 5-functions in Eq. (|A.8| ) 
and is expressed in terms of a single fiux line cumulant expansion of pi and fj^, in Fourier 
space, 

(l+m \ /l+m \ 

X (j){-ki, -qi) ■ ■ -(pi-ki, -qi)(f)rai{-ki+i, -qi+i) ■ ■ ■ (j)ra^{-ki+rn, -qi+m) , 

(A.ll) 

where the sum over is implied. Fal^^a^ is a connected correlation function of a single 
line density fluctuations pi and t[ , and in Fourier space is. 



(l+m \ /l+m \ 
Emm E^O 



A{piki,qi) ■ ■■P{ki,qi)ra,{ki+i,qi+i) ■ ■ ■ T^^{ki+m,qi+m))o ■ (A.12) 
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As in dynamics we have inserted an area factor A for convenience and dropped the single 
hne label 1. The vertex functions Ta^^^a^ can be systematically computed for any /,m. 
We calculate some of the vertex functions for both statics and dynamics in Appendices B 
and C. 

Within the Gaussian approximation we truncate r[(^, (^t] at quadratic order in the 
densities. We observe that T^^^^\k,q) = S°{k,q), T^^^'^\k,q) = T°^i,{k,q) and Ti^'^\k,q) 
are single line number and tangent static structure functions (with notation S°{k,q) and 
T"^fj{k, q) from Ref.|jl[] ), derived in Appendix B and discussed in the main text, 

rir'A«) = ^P:., (A.136) 
r<i-il(*,<,) = 0. (A.13c) 

The fact that T^a'^\kjq) and all the cross correlation functions vanish for the transverse 
part of the tangent density field rj is a result of the independence of rj from p. Math- 
ematically it is true to all orders because any cross correlation function of is a tensor 
with atleast one index a which must be carried by ka- This means that all the cross cor- 
relation functions are purely longitudinal and when contracted with (for every rJ) 
automatically vanish. 

Integrating over the auxiliary fields (p{k, q) and (prik, q) we obtain an effective Hamil- 
tonian expressed solely in terms of the physical number and transverse tangent density 
fields, 

n[^J] = i [^^(^,?)A-\^,Q)^(-^,-Q) + ^cA9RV(^.?)^n-^.-?) 

Jk,q ^ 

+^(k,q)U{-k,-q)^ , (A.14) 

where A{k,q) and Bab{k,q) are the static interacting but disorder-free 2pt-correlation 
functions of p and rJ, 

A(k,q)= ^oksTe/e 

q^ + iqBik)/kBTr 

BaBiKq)=PL^^. (A.156) 
e 
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and qsik) is the Bogoliubov spectrum given by Eg. ( [7.12]) in the main text. 

We use above Hamiltonian to derive the interacting static functions in the pres- 



ence of disorder, S'*(/c, = {p{k, q)p{-k, -q)) ^ and T^^^{k,q) = {Ta{k,q)n{-k,-q)) ^ = 
T'^ {k, q)Pab + "^sik, q)Pab- The longitudinal part of T^^ can be obtained from S^{k, q) by 
using flux line continuity Eq.( |A.^ ) . Averaging over the annealed density fields, followed 
by the average over the quenched disorder we obtain, 

S(k q)= ""o^BTkye ^ ^ ^ f Uoe/e \ ' 

q^ + {qB{k)/kBTY ' '\q^ + {qB{k)/kBTY) ' 

Tj(fc,5) = ^^, (A.166) 
T^{k,q) = ^Ss{k,q) . (A.16c) 



Eq.( |A.16a[ ) is in complete agreement with static limit of the interacting dynamic structure 



function obtained in the main text, Eq. (|7.11"^ . For short-range disorder and flux line 
interaction, F{k,q) = Aq and V{k) = (Pq/Stt'^, respectively above equations reduce to the 
results of Ref.BSl obtained using the boson mapping. 
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Appendix B. Calculation of the Nonlinearities in the Static Cumulant Expan- 



sion 



In this appendix we calculate the coefficients of the nonlinear interactions which appear 
in the cumulant expansion of the Hamiltonian in Appendix A. 

We begin with static problem and calculate the static single line correlation functions 
of the density fields p and rj, defined by Eq.( |A.12"D . It is convenient to work in the space 



of (/c, z). We evaluate the following static correlation function, 

V^aA^iki, Zi--- ki+^, Zl+rn) "4"*^^^^ ( 5Z ^ 

\i=l 

= i'^lA^^^ " "^a^^m, Zm)pikm+1, Zm+l) ■ ■ ■ Pikm+U Zm+l))o , (B.la) 

where the averages are performed with a Boltzmann weight of a single line and the super- 
script T extracts the transverse part of the average. 

Above average can be easily computed by using the following Lemma that holds for 
Gaussian averages, 



(Ai---A;e^)„ 



I 



e 



where, Xi and Y are Gaussian random variables with respect to the measure used. Above 
equation can be easily proved by rewriting the left-hand-side as an multiderivative with 
respect to / parameters of an exponential generating function, 

(Xi ■ --Xie^), = d^, ■■■d^, |{^^}=o (e^+'^^^Oo , (B.3a) 
= ■ ■ ■ d,, |{^.}=o e<^^>°/2+<^'^»^^>°+«'^^^^)^>°/2 , (B.36) 

where we used a property of averages over Gaussian variables, (exp(^))o = exp((^^)o/2). 
Performing the differentiation and evaluating the result at {pi} = we obtain Eq. ([B.2| ) . 
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Returning to the original problem, Eq.(B.l ) , we make the identification 

dziTaizi) and Y = -i Yl^^T ^» ' ^(^0^ obtaining, 



{Y')o = -{^\^}^^h-r{z,)j j , (B.4a) 

l-\-m Z+m 



l+m \"\ 
2 

''l+m \ Z+m 

^ fc, J C(0) - 2 ^ ■ - Zj) - C(0)] , (B.4c) 

where the r{zi) averages are isotropic, 

^{f{zi) ■ f{zj))o = C{z^ - Zj) , (B.5a) 



^iq{zi-Zj) 



(B.56) 



Using this expression we find. 



iq{zi-Zj) _ 1 

C{z, - z^) ~ C(0) = / , (B.6a) 



I '^j I 



(B.66) 



2e 

From Eq.( B.56 ) we observe that C{0) ^ L ^ A oo and therefore this first factor leads 
to a 5-function a,s L, A ^ oo and imposes momentum conservation in the xy-plane for all 
the correlation functions. 

The {XiY) and {XiXj) averages can be similarly evaluated, 

Z+m 

{X,Y)o = -id,^ ^(r(^^)^j • r{z,))o , (B.7a) 

Z+m 

= -z ^ /cj-5f(zi - ^j ) , (B.76) 



where, 



f (^i - = {r{zi) ■ r{zj))o , (B.8a) 
2 /"°° (ig sing(2;i — Zj) 



(B.86) 

- — sgn{zi - Zj) , (B.8c) 
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with sgn{z) = 1 for 2; > and sgn{z) = 



-1 for z <0. 



{XiXj)o = d,,d,.{r{zi) ■ r{zj))o , 

Jq 

= 25{zi - Zj) . 



(B.9a) 
(B.96) 
(B.9c) 



Combining above equations inside Eq.(B.l ) we find for odd m the cumulants vanish 
and for even m we obtain, 



{ki.Zi ■ ■■ki+m,Zi+m) 

e 



+A11 pair combinations] x exp 



kBT 
2e 



l+m 



^ ^ ki ■ kj I Zi 



i>j 



(B.IO) 



For Gaussian approximation that we are concerned with in the main text we obtain, 



r(i'i)ft^i,-^,^2) = o. 



(B.lla) 
(B.116) 
(B.llc) 



These expressions when Fourier transformed to q space give Eqs.(A.13 ) of Appendix A. 
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Appendix C. Noninter acting Flux-Line Liquid Dynamics 

In this appendix we apply methods similar to Appendix B to calculate the dynamic 
cumulants 

r2^(a^l,X2) = {pa{Xi)pb{x2))o , (C.l) 

where x = (/c, z,t) and a, 6 = 1, 2. The procedure is similar to the one used for statics. 
The averages are performed with the exponential of the single flux line dynamic functional 
Jo[f,r\, Eqs.( |3.9|) -(3.10 ), instead of with the static Boltzmann factor. As explained in 
the main text in the derivation of hydrodynamics we traced over the tangent density 
field and its corresponding response field. The longitudinal part of the dynamic tangent 
correlation function can be extracted from those of the number density field by using 
flux line continuity Eq. (|A.3| ) , while as for the statics, the transverse part of the tangent 
correlation function is probably not affected by the flux line interaction or disorder. 
In order to calculate the averages in Eq. (|C.1| ) we introduce "Rouse modes" , 

I 

r{q,t) = — I dzcos{pqz)f{z,t) , (C.2a) 
Jo 

r(z, t) = fo{t) + 2 J2 cos{pqz)r{q, t) , (C.26) 

where Pq = ^ and L is the length of the flux line. ro{t) denotes the center of mass 
mode. From the single line dynamic functional j7o[r, r], Eqs.( p.9D -(3.10 ), one can derive 
the following correlation functions for the "Rouse modes" , 

{r-{q',Oy{q,t)), = ^f^e-''^P>\dqq,dai,, (C.3a) 
2Lep^ 

(Kit) - (rlit) - r^(0))). = ^^^\t\6at , (C.36) 

0)r\q, = -±.6{t)e-^^P>\ 6qq,5ab , (C.3c) 

(^:(0)r^(t))o = -^e{t)dai> , (CM) 

{r:{0)rlit))o = . (C.3e) 
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Using the above equations we deduce for the segment-correlation function of the internal 
modes 



(K(z,t)-r"(z',0)] [r\z,t)-r\z',Q)\), 



^ . 2DkBT 2L ^ keT 
OabS ; 1^1 + — 2^ — ^ 

2L ^ knT 



^ (^cos{2pgz) + cos(2pqz')^ + cos{pq{z + z'))e~^^P^^^^^ 



(C.4) 



This has to be inserted in the corresponding expression for the single line correlation 
function S°{k, z, z' , t) = (^5 ^? ^'5 ^) 

rg)(fc, z', t) = {p{-k, z, t)p{k, z\ 0))o = exp |-lA;2([f(z, t) - r(;2', 0)]')„| . (C.5) 

Since we are dealing with a finite size system there is no translational invariance with 
respect to the 2;-axis, and therefore correlation functions like Eq. (|C.5|) depend explicitely 
on z and z' . But, except for the small region where z and/or z' are close to one of 
the ends of a flux line the second term in Eq.( |C.4| ) can be neglected. This is due to the 
rapid oscillations of cosine terms like cos(2]5q2). In the rest of this appendix we will neglect 
those contributions and hence get expressions which depend only on the relative coordinate 
{z-z'). 

For t > tRouse = ^ the summation ^^yi is rapidly converging. If t ^ ^Rouse the 
sums can be replaced by integrals 



where 



/(.) ^ i r ^ 



1 - e 



-yx 



COS X 



(C.6) 



(C.7) 



TV Jq 

The function f{y) can be written in terms of the incomplete Gamma-function r(a, x) as 



f(y) 



1 1 



2 4v^ V 2'% 
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(C. 



with the asymptotic representations 



f{y) 



for y > 1 , 
for y ^ 1 . 



(C.9) 



With the latter definitions we find for the two-point function 



r55(fc,^,t) = exp 



DkBTk' 



\t\ - k 



2/CbT 



Deltl 



(C.IO) 



Since /(O) = 1/2 we recover the static result of Appendix B as we must for the equal time, 



ti — tj = 0, correlation function. 



Next we consider the single line response function S^^\k, z, z' ,t) = F^^^ {k, z, z' ,t) = 
{p{—k, z,t)p{k, z',0))o, again in a regime where z and z' are not too close to one of the 
ends of the flux line. We find 



T^^^ik, z, t) = iDkBTe^{r{0, 0) ■ f{z, 

DkBTk"^ 



1 



exp 



t\~e^\z\f^ 
e V 



e{t)k'^DkBT 



1 



1 

L ^ M AnDeltl 



-z^/i4De\t\) 



DkBTk\ ,2^sT, .,rDe\t\ 
X exp <^ \t\ — k \z\f 



For the correlation function F]^]^ (/c, z, t) one gets 



(C.lla) 
(C.116) 



T^^^{k,z,t) = - {DkBTkf exp 



DkBTk'^ ,2kBT , , „ fDe\t\ 



t - k' 



\z\f 



^ {fiz, t) ■ f{0, 0))o - ^ ( {fiz. t) ■ r{z, t))o - (r(2, t) ■ f(0, 0))„ 



X (^(0, 0) ■ f{z, t))o - (r(0, 0) ■ f(0, 0))o 



. (C.12) 



All the terms inside the square brackets in Eq.( |Cl'2 ) vanish. (r{z^ t) •r(0, 0))o = as can be 
seen from Eq. (|3.10^ , and the equal-time correlators {f{z, t) ■ r{z, t))o and (r(0, 0) ■ r(0, 0))o 
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vanish by the causahty together with our choice of causal time discretization, Eq.( p.l4 ) . 
FinaUy, the remaining term {r{z, t)-f{0, 0))o(r(0, 0)-f'{z, t))o = 0, because it is proportional 
to 9{t)9{—t) which vanishes by the definition of 9{t) and by the choice of causal time 
discretization. Therefore, 

rfi^(A?,2,t) = . (C.13) 

We now observe that {k, z,t) and {k, z,t) satisfy the standard Fluctuation 
Dissipation Theorem (FDT), 



dtT^Sik, t) = ^'12 ik, t) - ^'12 ik, -t) 



(C.14) 



This is the reason why in the main text we called pi = p the response field of the physical 

(2) (2) 

density field P2 = p and called = and T21 = S" the correlation and response 
functions, respectively. 

Now we discuss the correlation functions in various limiting cases. For ksTk'^L / e <ti 
1, where the wavelength is much larger than the static transverse "radius of gyration" 



\plk^TT7[t of the fiux line, the dominant term in Eq. (|C.10|) is the first term. This term 
describes the center of mass motion of a single flux line. We flnd exponentially decaying 
correlation functions 



S°{k,z,t) UBTfc2L/e«i ^exp 



DksTk^' 



\t\ 



S {k,z,t) \kBTk^L/e<^i ~ 9{t) exp 



DkBTk"' 



(C.15a) 
(C.156) 



All the other terms are of order 0{k'^L) and can therefore be neglected. 

For ksTk'^L / e ^ 1, where the wavelength is much smaller than the "radius of gyra- 
tion", the internal modes, i.e. the second term in Eq. (|C.10|) , dominate provided t is not 
too large. For times t > tRouse the center of mass motion becomes the dominant contri- 
bution again, since the contribution from the internal modes scales as ^/i for large times 
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(note /(y > 1) ~ y/y). Hence we get for t < ^Rouse 

S%k, z, t) \k,Tk^L/e»i^ exp ^ \z\f (^^^ I (C.16a) 



S°{k, z, t) |feBTfc2L/e»i~ 0{t)k^ W — exp 



^2 



ADet 



e^^{-e^-^\z\f{^)\ (C.166) 



2^ 



Finally, we study the Fourier-transforms of these noninteracting structure and re- 
sponse functions to be used in the main text to calculate the interacting counterparts. 
In the center of mass limit one gets a Lorentzian shape for the correlation function 



(/(Q = 0) = /,^ci^/(^)), 



with the linewidth 



2LT (k) 

S"{k,q^O,u) \ksTk^L/e^l ~ ^ (C.176) 



rUk) = ^k'. (C.18) 



In the limit where the internal modes dominate we find 



4 /'°° I- r 

S°{k,q,uj)= ^ die'^* da;cos(^Vte)e-^'/^-^*"»(^) , (C.19) 

ksT sjT:k^ Jq Jq 

where we have introduced the scaling variables Cj = oj/V^, q = 2qe/ {ksTk'^), and = 
{kBT)'^Dk'^ /{Ae), the noninteracting dynamic relaxation rate. We have also defined a 
function m{x) by 

1 / 1 

m{x) =x + x—^T {--,— ] . (C.20) 



20F V 2' 4 

Expanding in powers of a) and we obtain 



e 4 

n=0 
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with the coefficients 



6„(g,,) = (,„«)= 2"+' (2n-H)!! f fe-V^ ^"- 1(2" + 2) ""^t^" (c.22) 
with the expansions 



where the coefficients are given by 



/"OO 

6„ = 2'^+i(2n + l)!! / dxe-^'/4m(x)-2'^-2 , (C.24a) 



n! 



The expression for S°{k,q^uj) can now be easily obtained from the above expression 
for S"{k^q,uj) and the FDT, Eg . (10.141) . In the limit of small frequencies one gets to 
leading order 

S°(i,9,a;)«^-|p[6o(«) + »i(<j)'i-f>2(<j)'i"] . (C.25a) 



^ 63(5) 

In the limit [uj 0) the Fourier transform S°{k,q,uj) can also be written in the 
Lorentzian form 

S^{k, q, u) ^ Mk, q) , (C.26) 

where we have introduced the quantities 

and 

T,,^,Uk,q)=TJ^-^^ (C.28) 
V h{k,q) 
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In the limit, 2qe/{kBTk^) <^ 1, As{k,q) reduces to 



Asik,q) = 2 



4 ^i-6(i^)^i 



keT 0F/c2 

k'^keT 
0re 



bi/bs 



3b[ 

261 2b3 



^1. 



(C.29) 



where 6„, b'^ are constants of order 0(1). The formula for Ag{k, q) is quite similar in form 
to the static structure factor of non-interacting lines 



In the same limit the linewidth reads 



q^ + {k^kBT/2ey 



(C.30) 



rsingie(A;,9) = r^ ( 



26i 263 



(C.31) 



Finally, we discuss the Fourier transform of the single line correlation function in the 
limit q = 

S^{k,q = 0,u)^ dz dte'''*S^{k,z,t) . (C.32) 
Jo J -00 

Upon introducing the scaling variables Cj = a;/(r^(l + k~^)), k^ = k'^LksT / {Ae), and 
K = r^^Rouse) where the latter variable gives the ratio of the Rouse time to the decay time 
of density fluctuations we obtain 



S{k,q = 0,uj) 



r,-(i + /c-2) 



F(k,(jj, k) 



(C.33) 



with 



F{k,u,K) = 



(C.34) 



c/r cos (tur) exp 



k- 



dyexp 



T 



( 



1 + /C-2 



rn 



1 + 
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The above scaling analysis of the single line correlation function shows that there are 
two time scales. First, there is the time scale for the dynamics of the internal modes 
given by tmtemai = l/r^. Then there is a time scale tRouse = -^^/(-De), which gives the 
crossover time above which the dynamics starts to be dominated by the center of mass 



mode. Furthermore, we have a length scale R(P = y '^^bTL describes the projected 

2D radius of gyration. 

For k <^1. i.e. for wave length much larger than the 2D projected Radius of gyration, 
the dynamics is determined by the center of mass motion, and we find for the scaling 
function 



F{k, Co, k) = ^^-^ , (C.35) 



l-e"^ 2 

which corresponds to Eq. (|C.17a[) , provided k is small, i.e. the Rouse time is small compared 
to the typical relaxation time for internal modes. 

For k ^ 1 the internal modes dominate the dynamic structure factor. In fig. ^ 
the scaling function F is shown for a fixed value of k = 1.0 and a series for k~^ {— 
0.0001, 0.5, 1.0, 2.0). For any value of k, and k the curves are rather well represented by 

^ (C.36) 

where the coeffiecients ao, ai and the exponent 7 depend on k and k. In the limit of small k 
the scaling function F turns into a Lorentzian with 7 = 2 and oq = 1.264, ai = 1.0, which 
just corresponds to the Lorentzian shape in the center of mass limit (compare Eq. (|C.17a| ) 
and Eq.( |C35| )). The exponent 7 decreases with increasing k. The best nonlinear fit for 
k~^ = 0.0001 gives ao ~ 2.93, ai ~ 5.93, and 7 ~ 1.344. Hence, the lineshape crossover 
can essentially be characterized in terms of the effective exponent 7. The typical linewidth 
r(/c, L) of the correlation function is given by the condition Cj = 1, i.e. 
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which describes the crossover from a dynamics determined by the internal modes to a 
dynamics governed by the center of mass motion. 

We note that S°{k, z, t) and S°{k, z, t) from Eqs. (|C.10|) , (C.ll )satisfy the static hmit 
sum rules, as they must, 

S°^{k,q) = S"ilq,t = 0), 
du) 



S'^{k,q,u), (C.38a) 
S"^{k,q) = S"{k,q,u = Q) . (C.386) 
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Figure Captions 

Fig. 1. (a) Mean-field phase diagram of type II superconductors, (b) Schematic picture 
of a phase diagram of high-Tc superconductors which includes effects of thermal 
fiuctuations and disorder. 

Fig. 2. Schematic picture of flux lines in high-Tc superconductor in the presence of pin- 
ning disorder (indicated by black circles). Conformation of the z-th line is de- 
scribed by a two-dimensional vector fi{z,t). 

Fig. 3. 2D-projected configuration of a diffusing fiux line for (a) t » tuouse and (b) 
t ^ t Rouse- In (ct) the progression in time is illustrated, with the inset showing a 
snap-shot configuration, "diffusion" in 2;-direction. For t ^ tuouse the flux line 
diffusion is dominated by the center-of-mass mode and the dynamics is that of a 
rigid rod. The dashed line shows the trajectory of the center of mass. In (b) the 
fliix line appears frozen with dynamic fluctuations (due to the center of mass and 
internal modes) small relative to its transverse size. The inset show evolution in 
time. 

Fig. 4. Scaling function f{x) describing single flux line diffusion due to internal modes, 
displayed on a double logarithmic scale. 

Fig. 5. Illustration of the interaction length Lj{k) deflned in the text. 

Fig. 6. Scaling function F for a fixed value of k = 1.0 and a series for (= 
0.0001 (solid), 0.5 (dotted), 1.0 (dashed), 2.0 (dot - dashed)). 
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